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ABSTRACT 

SZ clusters surveys like Planck, the South Pole Telescope, and the Atacama Cosmology Telescope, 
will soon be publishing several hundred SZ-selected systems. The key ingredient required to transport 
the mass calibration from current X-ray selected cluster samples to these SZ systems is the Ysz-Yx 
scaling relation. We constrain the amplitude, slope, and scatter of the Ysz-Yx scaling relation using 
SZ data from Planck and X-ray data from Chandra. We find a best fit amplitude of ln(_D^l52 / CYx) — 
-0.202 ± 0.024 at the pivot point CYx = 8 x 10"^ Mpc^. This corresponds to a Ysz /Yx-rdXio of 
0.82 ± 0.024, in good agreement with X-ray expectations after including the effects of gas clumping. 
The slope of the relation is a = 0.916 ± 0.032, consistent with unity at « 2.3(t. We are unable to 
detect intrinsic scatter, and find no evidence that the scaling relation depends on cluster dynamical 
state. 

Subject headings: clusters 
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1. INTRODUCTION 

The Sunyaev-Zeldovich effect — i.e. the predic- 
tion that the hot intra-cluster gas of massive galaxy 
clusters will scatter CMB photons to produce a 
unique frequency signature — traces back to 1972 
(jSunvaev fc Zeldovich I1972D. Despite early detections 
of this effect (|Gun fc Northoved 119761 ). it is only re- 
cently that instrumentation advances have made large 
scale SZ searches feasible. There are three such ongo- 
ing surveys: one carried out with the South Pole Tele- 
scope (SPT iCarlstromet al.l 1201 1[ ) one with the Ata - 
cama Cosmology Telescope fAC T iFowler et al.l 120071) 
and on e using the Planck satellite (jPlanck Coll aboration! 
I2011ef ). All three surveys have published initial chistcr 
samples (IVanderlinde et al.l 120101 : IMarriage et al.l 120111 : 
iPlanck Collaborationll2011d[) . and are expected to even- 
tually publish hundreds of SZ selected systems out to z = 
1 and beyond. If the scaling relation between the inte- 
grated SZ signal and cluster mass can be adequately cal- 
ibrated, then these cluster samples may be used to con- 
strain the growth of large scale structure in the universe, 
which can in turn provide critical con straints on dark 
energy models (e.g. iHolder et al.ll200il) . Indeed, both 
SPT and ACT have already derived cosmological con- 
straints from their current samples ("Va nderlinde et aP 
[2010: Sehgal et al. 2011: Benson et al. 2011^^ 

Unfortunately, at this time no single method provides 
masses of individual clusters with a required ^ 5% un- 
certainty. X-ray hydrostatic measurements are obser- 
vationally expensive and limited to a subset of relaxed 
systems. Numerical simulations suggest that hydrostatic 
mass est imates can be biased lo w anywhere between 10% 
to 30% (iNagai et al.. 2007: Piffaretti fc Valdarnini, 20081: 
iJeltema et al.l 120081: iLau et al.l 120091: iMeneghetti et al.l 
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120101 : iRasia et"al] I2012D . Weak lensing estimates are 
in principle unbiased, but specific implementations can 
result in w 5% — 10% biases depending on the fitting 
procedure (t ypically biased high when using NEW fit- 
ting formulaelBecker fc Kravtsovl2010HMeneghetti et al.l 
l2010HRasia et al.!l2012D .' and exhibit large intrinsic scat- 
ter. In addition, without the help from X-ray data, 
the derived SZ signal from the current generation sur- 
veys has a comparatively large scatter at a fixe d mass 
(jBenson et al.l 120 1 1! : [pTanck C oUabor ationll2011d( ). Mass 
calib ration using velocit y dispersion s is a lso possible 
(e.g. lEvrard et all l2008t iSifon et al.l |20T1) . but sub- 
ject to velocity bias. Therefore, the precise calibra- 
tion of the cluster masses from SZ surveys will be pos- 
sible only through a combined analysis of the SZ, X- 
ray, and optical data. Consequently, all three of the 
on-going SZ experiments are devoting significant effort 
for extensive X-ray and weak lensing follow-up of the 
SZ-se l ected systems (lAndersson et al.l 120101: iFolev et al " 
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IPlanck Collaboration! I2011bl |j: iMenanteau et al 
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Another possibility, however, is simply to transfer the 
existing calibratio n of the mass-proxy relations in the X- 
ray samples (e g. lArnaud et al.! 120071: lAllen et al.! 120081: 



Vikhlinin et~aIlT2009l: lArnaud et aLM2010l: iMantz et alj 



using the correlation of the cluster SZ parame - 



ters with the X-ray proxies (e.g. IBenson et al.! [20Tll) . 
An important step in this process is to verify a tight 
correlation between the integrated SZ s ignal, Ysz, and 
its X -ray analog , Yx = Tx x Mgas x ()Kravtsov et al.! 
l2006l see section 11.11 for more precise definitions of Ysz 
and Yx). In addition to transferring the mass cali- 
bration, the Ysz — Yx relation can be used to verify 
some key predictions of the structure formation theory 
which is used for applying the cosmological models to the 
cluster data. First, numerical simulations consistently 
show that both Ysz and Yx have low scatter at a fixed 
mass (» 10%o - 15%. INagail 120061: iKravtsov et al! 120061: 
iStanek et alllMol : iFabian et al.!l201l!: iKav et al.! 
Thus, we can expect the Ysz-Yx relation to have very 
low scatter, especially since these parameters approxi- 
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mate the same physical quantity, the thermal energy in 
the cluster gas. Second, these same simulations predict 
nearly self-similar evolution in the Y — M relation for 
both Yx and Ysz, so we expect the Ysz-Yx relation to 
have nearly the same normalization at both low and high 
redshifts. 

The Ysz-Yx r elation has been r eport ed in several re- 
cent works — bv lAndersson et aTl (120101) at high-z using 
SPT an d Chandra, and at low-z bv lPlanck Collaboration! 
(|2011af l using Planck and XMM-Newton. The results are 
generally in line with the theoretical expectations — Ysz 
follows Yx with a small overall offset. However, there are 
several issues which warrant further study. The overall 
offset between Yx and Ysz is expected because real clus- 
ters are non-isothermal and slightly non-uniform. Given 
the overall weighting of X-ray emission versus the p 
weighting of the SZ-signal, one naturally expects the two 
quantities to differ (see below for further details). Based 
on the analysis of XMM-Newton tempe rature profiles 
in a " representative" sample of clusters, lArnaud et al.l 
(|20T0t ) predicted Ysz/Yx = 0.924. The mean offset in 
the SPT /Chandra sample of lAndersson et all (|2010( ) is 
somewhat below this prediction, Ysz/Yx = 0.82 ± 0.07. 
The mean offset in the Planck / X MM-Newton sample is 
1 a ab ove, Ysz/Yx = 0.95 ± 0.03 (jPlanck Collaboration! 
I2011al) . While the difference in the Ysz/Yx ratio is rel- 
atively large, it is not signific a nt be cause of the large 
errors m the lAndersson et al.l ()2010l ) measurement, so 
additional measurements that reduce the statistical un- 
certainty of Chandra estimates of the Ysz/Yx ratio 
are of great interest. This is particularly true in light 
of known differences in instrument calibration between 
Chandra and XMM so having a low-redshift calibration 
of this ratio using Chandra is of significant interest. Fi- 
nally, the scatter in the Planck / XMM Ysz-Yx relation 
is 23% ± 2%, which is significantly larger than the ex- 
pected w 10% scatter, so confirmation of this observation 
is important. 

The goal of this paper is to measure the local Ysz — Yx 
relation in dependently using over lapping clusters in the 
Planck and lVikhlinin et all (j2009( ) Chandra sample. This 
is the first step towards C/iandra-based calibration of the 
local Ysz — M relation. We defer the actual calibration 
of the Ysz^M relation to an upcoming paper, where we 
simultaneously consider several the Ysz-M scaling rela - 
tion as calibrated by several groups ( Rozo et al.ll201l| ). 
In this context, we emphasize that the Chandra X-ray re- 
sults quoted here are derived from the gas temperature 
and gas density models for individual clusters taken from 
iVikhlinin et al.l(l2009D. C onseauentlv. this work and that 
of IVikhlinin et al.l (|2009f ) are fully self-consistent in all 
definitions. 

Throughout this work, we assume a flat ACDM cos- 
mology with flr n = 0.3 and h = 0.72 . This convention 
matches that of I Vikhlinin et al.l ()2009D . 



1.1. What We Mean by the Ysz/Yx Ratio 

The Ysz signal of a galaxy cluster is related to the 
total pressure of the intra-cluster gas within the cluster's 
volume via 



where Da is the angular diameter distance, and C is 
a constant that switches the units back and forth from 
M0keV to Mpc^. The constant C is given by 
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where ctt is the Thompson cross-section, and Pgas/rie = 
1.149 mp for fully ionized plasma with the cosmic He 
abu ndance and abundan ce of heavier elements set at so- 
lar. lArnaud et al] (I2010D assume a heavy-element abun- 
dance that is 0.3 solar, which differs from the above value 
by less than 1%, a difference that is completely negligible 
for our purposes. 

From equation [TJ it is apparent that the combination 
D\Ysz will scale with CYx where Yx = Mgas x Tx 
(see section 13.3! for our precise definition of Mgas and 
Tx)- Thus, we expect that the Ysz/Yx ratio is a di- 
mensionless number that is close to unity. Throughout 
this work, we will often speak about the Ysz/Yx ratio 
in this fashion. In practice, however, Ysz and Yx are 
cluster observables that are related via a probability dis- 
tribution P{Ysz\Yx)- Here, we assume P{Ysz\Yx) is a 
log-normal distribution, which takes the form 

(3) 



(In {D\Ysz) \Yx)^a + a\nCYx. 



When we speak of "the Ysz/Yx ratio", we mean the 
quantity (\xi[D\Ysz) \ Yx) — In CYx evaluated at the 
pivot point of the sample. Note that if we select units 
such that CYx = 1 at the pivot point, then the Ysz/Yx 
ratio is just the amplitude a. We will work primarily 
in units of CYx = 8 x 10~^ Mpc^, corresponding to 
the pivot point for the Chandra- Planck cluster sample 
(see section [3]) . For brevity, we will often simply write 
\n{Ysz/Yx) — as opposed to the fuU {in D\Ys z\Yx) ~ 
In CYx expression — to denote the Ysz/Yx ratio. 

2. X-RAY PREDICTIONS 

Equation [1] allows one to compute the Ysz signal of 
a galaxy cluster based on the gas and temperature pro- 
files of galaxy clusters. With sufficiently high quality 
X-ray data, both of these profiles can be observation- 
ally constrained, allowing one to directly predict the 
Ysz signal of a g alaxy cluster from t he X-ray data. 
Using this method, lArnaud et al.l ()2010[ ) predicted that 
D\Ysz/CYx = 0.924 or \n{Ysz/Yx) = -0.08 when Tx 
is defined withi i i a [0 .15, 0.75]i?5oo aperture. Using the 
IVikhlinin et al.l ()2006D Chandra results for local relaxed 
clusters, we predict nearly the same ratio. Specifically, 
the mass-weighted temperature 
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D\Ysz 
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is related to the spectroscopic temperature Tx via 
Trnw/Tx = 0.928, where Tx is measured in the [0.15 - 
1] i?5oo aperture. Using the definition for Tmw, we see 
that we can rewrite equation [T] via 

D\Ysz = CMgasT„^^, (5) 

from which we find \n{Ysz/Yx) = —0.07. 

It is potentially a concern that the lArnaud et al.l (|2010[ ) 
work is based on a representative cluster sample of X-ray 



Ysz~Yx from Planck and Chandra 
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luminous clusters while IVikhlinin et alj ()2006D include 
only relaxed, cool core clusters. However, empirically, 
there is very little difference in Ysz/Yx for relaxed and 
unrelaxcd clust ers .iudging; either fr om the X-ray predic- 
tions (Fig. 9 in lArnaud et a l.' 2010) or the actual mea- 
surements (see Fig. 4 in fPlanck C oll aborationll20lTaL and 
our results below). Another point of concern is slightly 
different apertures for Tx used in the XMM and Chan- 
dra samples. One expects that mean (emission-weighted) 
temperatures in the [0.15 — 1] i?5oo and [0.15 — 0.75] i?5oo 
annuli should be close because only a small fraction of 
the total X-ray flux comes from outside 0.75i?5oo. In- 
deed, the rA-(f0.15, 0.75] j^5 00)/ry(f0.15, l]i ? 500) r atios 
in the Chandra sample of IVikhlinin et al.l (120091 ) are 
within ±0.015 of a mean value of 0.985. An analysis 
of the deeper pointings among the ESZ-XMM sample 
(G.Pratt, private comm.) shows a similar value, even 
tho ugh the temperatu re ratio for the REXCESS sam- 
ple (jPratt et al.ll2009[ ) was somewhat lower on average 
(« 0.95), and showed a larger spread 0.85 — 1. We pro- 
ceed accordingly, and ignore the difference in the defini- 
tion of TX between the two works. With this caveat in 
mind, we conclude that both XMM and Chandra predict 
ln(Xsz/yx) — —0.075 because of non-isothermal cluster 
temperature profiles. 

There is, however, an additional complication that 
modifies these expectations. Specifically, when one esti- 
mates Mgas from X-ray data, one uses the fact that the 
surface brightness S cx p^^^ to estimate (pgas): averaged 
in radial annuli. For a uniform, spherically symmetric 
distribution (p^as) — (Pgas)"^ , so one can recover the gas 
density profile. In practice, however, one will typically 

have (p^as) > {Pgas)'^ ■ The clumping factor is defined 
via 

Q' = (6) 

{Pgas} 

SO that Mf^l ^ QMllT- By definition, Yx = Mf^'^Tx, 
so the Ysz to Yx ratio takes the form 

DaYsz _ 1 Tmw 
CYx ^ Q Tx' ^ ^ 

In practice, the clumping factor Q depends on radius, 
and the factor 1/Q characterizes the mean clumping cor- 
rection over the radius R^oo- Based on this discussion, 
we modify our X-ray expectation by a factor of Q~^. If 
there is no substructure masking in the X-ray analysis, Q 
can be very large (e.g. Q — l.lG. lMathiesen et al.lll999l ). 
though analysis of numerical simulations that include 
substructure masking achieve significantly lowe r values 
of Q. Based on the analysis of lNagai et al.l (j2007| ). we set 
Q — 1.05 ± 0.05, where the error is set to the size of the 
correction to account for systematic uncertainties in the 
simulation physics. The corresponding expectation for 
the Ysz/Yx ratio is \ii{D\Ys z / CYx) = -0.125 ± 0.05. 

Finally, we note that Tx measured in the X-ray can 
be subject to a systematic uncertainty, related to, e.g., 
absolute calibration of the X-ray telescope effective area. 
Assuming for simplicity that this results in a uniform 

bias of observed temperatures, x^^^"^ — AT^'"°'', the 



measured ratio of the y's can be factorized as 

DaYsZ _ 1 T^nw /ON 

CYx ~ AQ Tx' ^ ' 

Discussion of the calibration-related uncertainties in Tx 
is beyond the scope of this paper and we proceed assum- 
ing A = 1. However, it is important to keep in mind a 
possibility of a few percent systematic error in the X-ray 
temperature measurements with the current calibration 
of the XMM-Newton and Chandra instruments. To give 
a handle of possible deviations of A from unity, we note 
that within the sample we consider in this work, the aver- 
age ratio of XMM and C/iandra-measured temperatures 
for the same cluster is \n{TxMM /Tchandra) = -0.08 (for 
a detailed comparison of XMM and Chandra data on a 
cluster-by-cluster basis see Rozo et al. 2012, in prepara- 
tion) 

3. DATA 

3.1. Cluster Sample 

The cluster sample employed in this work is a 
subset of the Early SZ (E SZ) sample published in 
iPlanck CoUaboratioiJ ()2011d[ ). The ESZ sample is com- 
prised of 189 cluster candidates identified in a blind, 
multi-frequency search of the all-sky maps obtained from 
the first ten months of observ ations of the Planck mission 
(jPlanck Collaborationl[2011en . The threshold for cluster 
detection was set at a signal-to-noise (S/N) > 6. This 
ESZ sample was searched for archival XMM data, and a 
sub-sample of 88 galaxy clusters with XMM data were 
identified. Of these, 62 clu sters were processed in time fo r 
publication, as detailed in IPlanck Collaboration! ()2011al ). 
Starting from this sub-sample of 62 galaxy clusters, we 
identified 28 clusters with Chandra data that had al- 
re ady been analyzed us ing the data reduction pipeline 
of IVikhlinin et all ()2009D . The analysis presented here is 
based on this sub-sample of 28 galaxy clusters, which is 
summarized in Table [H 20 of these are relaxed galaxy 
clusters, and 8 are merging systems. A plot of the Ysz 
vs Yx data for the galaxy clusters is shown in Figure [TJ 

3.2. Planck Data 

For each of our galaxy clusters, (jPlanck Collaboration! 
I2011al) estimated the integrated Compton parameter Ysz 
using a m atched-filter al g orithm . The matched filter as- 
sumes an lArnaud et all (|2010l ) pressure profile, which 
is used to predict the CMB temperature perturbation 
ST{9) of each galaxy cluster in each of the Planck fre- 
quency bands. In order to perform the fit, the clus- 
ters are first assigned a putative R500 radius based on 
th e XMM X-ray data and the M-Yx scaling relation 
of lArnaud et al.l ()2010[ ). The amplitude of the matched- 
filter, which can be expressed in terms of the integrated 
Compton parameter Ysz, is then fit using the temper- 
ature data within a projected aperture of 5i?5oo. We 
not e that the shape of the match ed filter is held fixed in 
the IPlanck Collaboration! ()2011a|) analysis, but their sys- 
tematics tests demonstrate that varying the shape of the 
filter within tolerable levels does not significantly impact 

^ The quoted ratio docs not correct for the difference between the 
temperature definitions used in the relevant XMM and Chandra 
works as noted earlier. 
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Fig. 1. — The relation between the Planck D^Ygz and Chandra 
CYx measured in the same aperture, -R500 provided by XMM- 
Newton mass estimates. Blue and red crosses repr esent dynamical 
relaxe d and merging clusters, respectively (per thc lVikhlinin et al.l 
II2009I ) classification). The solid line shows a 1-1 relation, while 
the dotted line is our best fit relation. The amplitude at the pivot 
point corresponds to D^^Ygz/CYx = 0.82. 

the recovered scaling relations. The median error in the 
recovered Ysz values for our cluster sample is 11%, and 
systematic uncertainties due to instrumental beam size 
and radio point source contamination are all expected 
be at the few percent level a t most. Note that because 
iPlanck Collaboration! (|2011aD assumed h = 0.7 while we 
as sume h = 0.72, we must res cale the D\Ysz values 
in IPlanck Collaboration! ( )2011al) appropriately. Further, 
the Planck team uses the redshifts listed for the clusters 
in the NED and SIMBAD databases while we use the red- 
shifts of the dominant central galaxies. We believe the 
latter approach is more robust because it allows one to 
avoid contamination by projection effects in some cases. 
Therefore, we rescaled all angular diameter distances to 
our adopted z values. The cluster redshifts, along with 
the rescaled Ysz values, are detailed in Table [T] 

3.3. Chandra Data 

Our analysis is based on the sub-sa r nple of ESZ galaxy 
clusters from IPlanck Collaboration! (|2011d[ ) that have 
Ch andra data tha t was reduced as part of the work 
in IVikhlinin et aTl (l2009f). described in full detail in 
iVikhlinin et al.l (120051 120061) . Here, we provide only a 
brief summary. We note that no t all galaxy clusters in 
our analysis are included in the IVikhlinin et al.l (j2009l ) 
sar nple. These are the clusters which were analyzed 
bv IVikhlinin et all ()2009[ ) but not included in the flux- 
limited sample used in their cosmological analysis. The 
full list of clusters employed in our analysis is presented 
in Table m 

We use the available Chandra data to estimate Yx = 
Mgas X Tx for each of the clusters in our sample, 
where Mgas is the total gas mass with i n the radius 
i?50o employed by IPlanck Collaboration! ()2011a!) . and 
Tx is the spectroscopic temperature within an annulus 
R e [015, l]i?5oo- To estimate Mgas, we rely on ana- assigned observational error. A similar expression holds 



lytic fits to the gas density profiles of individual clus- 
ters. The fits were based on the X-ray brightness profiles 
in the 0.7-2 keV energy band measured with Chandra 
and ROSAT PSPC (see IVikhlinin etlll (pOOl ) for de- 
tails). All detectable X-ray substructures were masked 
out from the profile measurements to minimize the im- 
pact of gas clumping. The X-ray fiux in the 0.7 — 2 keV 
band is very insensitive to the plasma temperature so 
long as T > 2 keV, as in the case of the sample un- 
der consideration. Nevertheless, we apply the appropri- 
ate temperature and metallicity corrections to the X- 
ray profile to cstim a,te M„as. The mean cluster tem- 
perature reported bv IVikhlinin et all ()2009D corresponds 
to the single-temperature fit to the X-ray spectrum in 
the 0.6 — 10 keV band extracted within an annulus 
[0.15, l]i?5oo- The median error in the recovered Yx val- 
ues for our cluster sample is 3%, nearly identical to the 
median error in the correspo nding XMM X-ray data in 
IPlanck CollaborationI (|2011aD . 

3.4. On the Importance of Apertures 

As already mentioned in section [1.11 we ignore any ra- 
dial dependence on the Ysz~Yx scaling relation. In other 
words, we are explicitly assuming that the parameters a 
and a are insensitive to the radius at which both Ysz 
and Yx are evaluated. This is clearly not true over arbi- 
trarily large radial ranges, but we expect any evolution 
over the range of radii probed to be at the < 1% level, 
so long as both Ysz and Yx are always measured within 
the same aperture. S i nce we rely on the SZ data from 
IPlanck CollaborationI (j2011aD . we have to use the same 
metric apertures adopted in that work. Specifically, to 
estimate Yx, we first measure Mgas within the assigned 
i?50o aperture, and multipl y this by the es ti mated tem- 
perature Tx as quoted in IVikhlinin et all ()2009[ ). We 
do not remeasure Tx since the X-ray temperature is ro- 
bust to modest changes in the assigned i?5oo radius of a 
cluster. A relation between P/ancfc-mcasured Ysz and 
the Chandra Yx values is shown in Fig. [TJ There is an 
obvious, tight correlation that is nearly linear, with no 
obvious offset between merging and relaxed systems. 

4. LIKELIHOOD MODEL 

Because wc will be assuming log-normal models 
throughout, we begin by defining and ysz via 

CYx 



Vx = hi 



Vs 



:ln 



8 X 10-5 Mpc^ 

D\Ysz 
8 X 10-5 Mpc^ 



(9) 
(10) 



We divide both Yx and by 8 x 10"^ Mpc^ so that 
both quantities are still measured in the same units, 
while simultaneously shifting the pivot point of the ysz- 
yx scaling relation to ensure the amplitude and slope of 
the best-fit Ysz-Yx relation are nearly uncorrelated. 

For each cluster, we assume that the observed y^ value 
is given by 



(11) 



where is the true value, and 5y is a random gaussian 
fluctuation of zero mean and = cr^, where ax is the 
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for Usz in terms of the true SZ signal ysz- The true val- 
ues ysz and are themselves related via a probability 
distribution P{ysz\yx) that characterizes the scaling re- 
lation. We assume this is a log-normal distribution of 
mean 



{ysz\Vx) ^a + ay^ 



(12) 



and variance ct^ i , so that ct<j,u is the intrinsic scatter 
in the Ysz-Yx relation. We assume this scaling relation 
is rcdshift independent, so a, a, and <7sz\x are all redshift 
independent. Note that because we have chosen our units 
to ensure that = at CYx = 8 x 10"^ Mpc^, the 
amplitude a is the Ysz/Yx ratio at this scale. 

The probability that a galaxy cluster has observed val- 
ues yx and y^z can be written as 

P{ysz,Vx)^ / dy^dysz P{ysz,Vx\ysz,yx)P{ysz,yx) (13) 



dyxdysz P{ysz\ysz)P{yx\yx)P{ysz\yx)P{yx)-,{'^^) 



where we arc using the fact that the observational errors 
in yx and ysz are uncorrelated. In principle, there can be 
a small amount of covariance induced by aperture effects. 
We discuss this effect in more detail below, but the main 
conclusion is that the effect is negligible for our purposes, 
and so we will ignore this effect in our discussion. 

Note the probability P{ysz,yx) depends on the distri- 
bution P{yx) of the cluster sample under consideration, 
which depends on cluster selection effects. We consider 
two possibilities, namely a flat prior for the value y^-, 
and a power- law distribution in Yx, corresponding to 
P{yx) = exp(— /3j/a;). The slope f3 can then be set to 
the expected slope for a WMAP7 cosmology. In prac- 
tice, neither is exact, but it can help give us a sense of 
the level of systematic uncertainty in our data due to 
selection effects. 

Setting P{yx) to a constant, integration of equation 
[HI allows us to conclude that P{ysz,yx) is a Gaussian 
distribution 



Piysz,yx)^Gi6;<jl,) 
1 



: CXp 



2aL 



where 



S^ysz - {a + ayx) 



2 2 

-a 



(15) 
(16) 



(17) 
(18) 



For a power-law distribution P(]jx) oc exp(— /Jy^,), we 
find 



P{ysz,yx) oc G{6;(T'}^t) 



X exp 



1^2 



{yxfsz + yxfx) + Ti^'^totfxf.i 

2 



(19) 



where yx is the yx value expected from y^, from naively 
inverting the scaling relation, i.e. 

yx = a-i(y,, -a). (20) 

/a; and jsz are the fraction of the total variance due to 



and ysz respectively, 

^ SZ — 



9 ? 



sz\x 



<^tot 



(21) 
(22) 



Note that for /3 = (i.e. a flat prior), we recover equa- 
tion I16[ as we should. For our fiducial analysis, we will 
assume /3 = 1.8, as appropriate for a logarithmic slope 
of the cumulative halo mass function N{M) oc for 

3 /5 

7 = 3, and a self-similar slope M oc Y^ . We have ver- 
ified that setting /3 ~ 0.0 does not impact our result at 
any significant level, which suggests our results are not 
sensitive to the details of the cluster selection function. 
This is in good agreemen t with the more detailed treat - 
ment of selection effects of lPlanck CoUaborationl (j2011aD , 
who find that selection effects only impacts the slope of 
the best-fit relation, and then only at the Aa = 0.01 
level, a systematic that is three times smaller than the 
statistical uncertainty. 

We assume that all galaxy clusters represent in- 
dependent realizations of the probability distribution 
P{ysz,yx), so that the likelihood of the entire data set 
is simply the product of the individual likelihoods. The 
posterior is 

C{a,a,(Jsz\x\ysz, yx) oc Po{a,a,crs^lx) P{ysz\yx) 

clusters 

(23) 

where Pq is the prior, and the product is over all clusters. 
Note that while we are using agz\x as a model parameter 
in order to facilitate physical interpretation of our results, 
it is the variance cr^^j^ that is more physical, so we adopt 

a fiat prior on cr^^|^, corresponding to 

Po(a, a, asz\x) « asz\x- (24) 

4.1. Aperture Correlations 

As detailed above, we have assumed that both Yx and 
Ysz are estimated within the same aperture, and we have 
explicitly assumed that the Ysz/Yx ratio is radius in- 
dependent. This guarantees that the mean relation is 
always adequately sampled. However, any noise in the 
estimated aperture i?5oo necessarily induces correlated 
noise in the estimated Ysz and Yx values for each clus- 
ters, and this covariance should in principle be included 
in the analysis. In practice, however, this is difficult 
to achieve because the Yx a perture adopted is that of 
iPlanck Collaboration! (|2011aD rather than the i?5oo aper- 
tu re that we would assign based on the Yx measurements 
of iVikhlinin et all (|2009n . This makes developing a self- 
consistent model that accounts for this covariance diffi- 
cult. We can demonstrate, however, that we expect this 
effect to be negligible. 

To sec this, we go back to equation [TT] Wc now let yx 
be the true yx value at the true radius -R500, while yx{R) 
is the true yx value at radius R. Assuming a power-law 
scaling 

yx{R) = yx+lx\n[-B—]. (25) 



Ri 



500 



and letting the observed i?5oo radius be a power-law as 
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a function of the observed Yx , one has 
In I I = T{yx - yx) 



R. 



500 



(26) 



where r is some proportionahty constant. Our modified 
equation [TT] is 

y. = yAR) + Sx, (27) 

so putting it all together, we arrive at 

yx = yx + — Sx. (28) 

1 - IxT 

The gas density profile of galaxy clusters is such that 
jx = 1-2, and for a self-similar scaling, R oc M^/^ oc 
Y^^"", so T = 0.2. Thus, if we self-consistently estimated 
the radius i?5oo from Yx , we would find that the variance 
in Yx increases by w 25%. Since the variance in Ysz 
is much larger than that in Yx, this is not a source of 
concern. 

What about the impact these fluctuations have on 
Igz? A similar argument as the one above results in 



Vsz = ysz 



IszT 



1 



Is 



-Sx + Ssz- 



(29) 



Thus, the fluctuations in y^^ are now explicitly corre- 
lated with those in y^. Note, however, that ^ Ssz 
because the SZ measurements are noisier, and that the 
fluctuations are themselves down- weighted by a factor 
« 0.3. This is a very minor perturbation to the variance 
in ysz, and the correlation coefficient between yx and ysz 
induced in this manner is very small indeed. To first 
order in JszT and ^xT we find 



IszT- 



0.06. 



(30) 



Thus, we expect the covariance induced from a self- 
consistent estimate of i?5oo to be negligible in this case. 
This independence is further strengthened by the fact 
that i?5oo is not estimated from the Chandra data itself, 
so while we cannot adequately model the effects of aper- 
ture choice on the variance of our measurements, we do 
not expect these to impact our results at any significant 
level. 

5. RESULTS 

We evaluate the posterior on a grid in the 3- 
dimensional space spanned by the parameters a, a, and 
<^sz\xy where the range of parameters is chosen to as 
to ensure that it does not impact our results. We use 
a e [-0.35,-^0.05], a e [0.8,1.1], and asz\x G [O,O.25]0 

Figure [3] shows our best fit scaling relation. The best 
fit amplitude (i.e. the Ysz/Yx ratio), slope, and scatter 
are 



a = -0.202 ±0.028 
a = 0.916 ±0.035 



z\x ' 



: 0.082 ±0.035. 



(31) 
(32) 
(33) 



^ In section [6l we r epeat o ur analysis using data published in 
IPlanck Collaboration! Il2011al) . Wiien doing so, we always make 
sure tile extent of our multidimensional grid completely covers the 
high-likeliliood regions in parameters space. 




Fig. 2. — The distribution of the log-likelihood value from 
our cluster sample. Wc use our best fit model to evaluate the 
log-likelihood of each clusters, and a Gaussian Kernel Density 
Estimator (KDE) to estimate the density distribution. The black 
line uses Silverman's rule of thumb to set the width of the KDE 
(i.e. the width of the Gaussian is set to l.SAa/N'^'^ where N is 
the number of samples). The dashed blue curve halves this width, 
and is obviously under-smoothed. There are no glaring outliers. 
The most discrepant cluster is A2163 (solid red line), which has a 
low likelihood because of it's extremely large mass. 



In all cases, the quoted values are the mean of the dis- 
tribution, and the error is the standard deviation of the 
marginalized distribution. The slope is marginally con- 
sistent with a = 1 (pa 2.3(7 offset). 

We have tested the sensitivity of our result to outliers 
by performing jackknife resamples of our galaxy clusters, 
and repeating our analysis. The resulting mean and stan- 
dard deviation are 



a =-0.202 ±0.024 
a = 0.916 ±0.032 
,L, = 0.083 ±0.021 



(34) 
(35) 
(36) 

The central values and error estimates are essentially 
identical to those derived from our full likelihood anal- 
ysis, demonstrating that there are no large outliers. 
This can be furthered evidence by evaluating the log- 
likelihood of each individual cluster. The distribution of 
log-likelihoods is shown in Figure^ The least likely clus- 
ter is A2163, reflecting the extremely large mass of this 
system. By the same token, A2163 signiflcantly increases 
the lever-arm with which the slope of the Ysz~Yx rela- 
tion can be constrained, and has an important impact on 
the flt. Removing A2163 from the flt results in a slope 
a = 0.906 ± 0.046. The error increases significantly, but 
the deviation of a from unity is still 2a. 

We test whether the scatter of the Ysz^Yx relation 
is consistent with zero or not by focusing not on the 
scatter, but on the variance, cr^^i^, which is the more 

physical parameter. The marginalized likelihood for the 
variance cr^^|^ is shown in Figure [4l We see that the 
data is consistent with zero variance. We can place an 
upper limit on the variance, cTj^i^ < 0.021 (95% CL). 
corresponding to (7sz\x l£ 0.144. 
We also explicitly consider the Ysz-Yx relation with 
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Fig. 3. — Best-fit parameters for the Yg^-Y" x scaling relation. Black cur ves are obtained using Chandra data and the cluster sample 
in Table [T] Red and blue curves arc using the IPlanck Collaboration! Il2011al) data. The red curves show the results when we restrict our 
analysis of the Planck Collaboration (2011a) data to the same galaxy clusters that constitute our Chandra sample. The blue curves include 
all galaxy clusters in Pl anck Co llaboration ( 2011al ). All confidence regions are 68% and 95% likelihood contours respectively. 



a fixed slope of unity (a = 1). The scatter is some- 
what increased to (Tgz\x ~ 0.120 ±0.033. The uncertainty 
in the amplitude also increases, reflecting the additional 
intrinsic scatter required to fit all the data. We find 
a = —0.201 ± 0.032. In this context, we also investigate 
whether the resulting amplitude is sensitive to the dy- 
namical state of the galaxy clusters. We split the clusters 
depending on whether they are r elaxed or unrelaxed, us - 
ing the same criteria employed in lVikhlinin et all ()2009D . 
The best fit amplitude for the relaxed and unrelaxed 
cases is a = -0.17 ± 0.04 and a = -0.26 ± 0.06, con- 
sistent with no systematic offset between the two. 

6. DISCUSSION 

6.1. Comparison to Previous Work 

We compare our best fit Ysz-Yx scaling relation 
to other empirical investigations of the Ysz/Yx ratio. 



There are several analysis of the SZ and X-ray prop- 
erties of galaxy clusters, b ut most cannot he stra ight 
forwardly compare to ours. iBonamente et al.l (|2008[ ) do 
not explicitly consider the Ysz-Yx relation, but sim- 
ply note that an isothermal model with no gas clump- 
ing gives a good fit to their data. This model corre- 
sponds to \'n{Ysz/Yx) = 0, which is in strong conflict 
with our results. We note, however, that they used sig- 
nificantly smaller apertures (i?250o)! ^-^d that their SZ 
and X-ray estimators are explicitly linked in the analy- 
sis: the SZ data helps constrain the X-ray data model 
and vice-versa, and the isothermal (3 model effectively 
acts as a prior relating the two data sets. Thus, it is 
difficult to fairly com pare their an a lysis t o ours. 

On a similar vein, IComis et al.] (I2OIII) used the AC- 
CEPT cluster catalog jfCavagnolo et al.l 120091) to probe 
the Ysz^Yx relation, but within an overdensity A = 
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a 




0.00 0.02 0.04 0.06 
Variance( InCD^^Ysz) I Y^) 



0.08 



Fig. 4. — Likelihood distribution for the intrinsic variance in 
lnYs2 at fixed Yx- Black, red, and blue curves are for the Chan- 
dra-Planck data, the XMM— Planck data restricted to the subset 
of clusters with Chandra data, and the full XMM-Planck data set. 



2500. In addition, iComis et al.l ()2011[ ) do not exclude 
the core when estimating Tx , which is expected to intro- 
duce scatter, and to systematicahy bias the Yx measure- 
ments relative to what would be measured using core- 
excluded temperatures. Indeed, IComis et al.l ()2011|) find 
both higher scatter and a higher Ysz/Yx ratio, as we 
would expect. A more quantitative comparison with our 
results , however, is difficult. 

The IPlanck Collaboration! (|2011aD analysis is the one 
that is closest in spirit to ours. Further, our analysis 
relies on their SZ-data. Switching the best fit ampli- 
tude reported in that work to the pivot point adopted 
here, their best fit scaling relation corresponds to a = 
-0.037 ± 0.023 and a = 0.96 ± 0.04. This amplitude is 
significantly higher than ours at more than 4.5(7. Their 
slope is steeper, but consistent with ours. There is 
also tension between the recovered scatter of their anal- 
ysis and ours, with the £Vz^_^= 0.228 value reported 
bv iPlanck CollaborationI ()2011al) being ruled out by our 
data at very high significanceQ 

What drives the difference b etween our results and 
those of IPlanck CollaborationI ()2011a|) ? Two obvious 
possibilities that we can readily explore are the fitting 
method used to recover the scaling relation, and the fact 
that we use a subsample (28 out of 62) of the galaxy 
clusters employed in the IPlanck CollaborationI (|2011a[ ) 
analysis. 

We test the importance of the fitting method by re- 
fitting the IPlanck CollaborationI ()2011aD data with our 
likelihood method. We find a = -0.032 ± 0.028, a = 
0.917 ± 0.039, and a^zlx = 0.182 ± 0.025. For this 
analysis, we have used the d ata exactly as reported in 
IPlanck CollaborationI ()2011a|) . that is, we use their red- 
shifts, and we assume h = 0.70, though we note that 
the h value impact the amplitude a only at the 1% 
level. The onl y exception is cluster A2034j for which 
the redsh ift in IPlanck CollaborationI (|2011a[ ) is clearly 
incorrect (jPlanck Collaborationll20iral quote z = 0.151, 



while we use the cD galaxy redshift, z — 0.113, a very 
significant difference). These results are illustrated with 
the blue lines in Figure [3l Our best fit amplitude is 
in agre ement with that quoted in IPlanck CollaborationI 
(|2011al) . though our best fit value for the intrinsic scat- 
ter is lower, with Aa^^ix ~ 0.04 (partly because we have 
correct the redshift of cluster A2034), and our best fit 
slopes are flatter. Despite the differences noted above, 
we can conclude that the fitting method is not respon- 
sible for th e difference between our results and those of 
IPlanck Collaboration (2011ai) . 

We test the impo r tance of using a subsample of 
Planck Collaboration (l2011af) data by analyzing the 
Planck CollaborationI (l2011aD data, but restricting our 



^ The number quoted in IPlanck CollaborationI BOllal) is 0.10 ± 
0.010, but they quote scatter in dex, while we quote the scatter in 
\nYsz. 



analysis to the clusters employed in our Chandra anal- 
ysis. In this case, we find a = —0.166 ± 0.035, a — 
0.912 ±0.043, and a,^\^ = 0.132 ±0.036. We see that the 
sample selection plays a significant impact on the results 
of the recovered Ysz/Yx ratio for the XMM data. The 
higher normalization of the Yx — Ysz relation from the 
full Planck sample must therefore be driven by the re- 
maining 34 clusters not in the Chandra subsample (most 
of which are at z > 0.1). Note that while we have 
not attempted to rep l icate t he Monte Carlo analysis of 
IPlanck CollaborationI (|2011a| ) to try to estimate selection 
effects, their own tests suggest that there are negligible, 
having no impact on the amplitude of the Ysz~Yx re- 
lation, and impacting the slope only at the Aa = 0.01 
level, much less than the statistical uncertainty in the 
measurement. 

Even when restricting ourselves to the same galaxy 
clusters, however, the Chandra data results in a lower 
Ysz/Yx ratio. While this difference is < la in terms of 
the statistical error bar, we note the statistical signifi- 
cance is significantly higher because both data sets share 
the same Ysz data, and it is the Ysz error that domi- 
nates the error budget. More specifically, the fact that 
both clusters use the exact same Ysz data implies that 
the the difference must be driven entirely by systematic 
differences in Yx between the two data sets. 

Aside from the differences in amplitude, it is interest- 
ing to note the difference in the intrinsic scatter esti- 
mate between the two data sets. The simplest interpre- 
tation of this result is that there a re unidentified sources 
of stochasticity in the IPlanck Col laboration (2011a|) X- 
ray analysis that result in an overestimate of the intrinsic 
scatter of the Ysz/Yx relation of galaxy clusters. 

Finally, it is worth noting that all of our fits consis- 
tently return a s lope a » 0.91, compared to th e a = 0.96 
slope reported in lPlanck CollaborationI (|2011a| ). This dif- 
ference appears to be due to the choice of fitting method: 
fitt ing our data using the same B CES method employed 
bv IPlanck CollaborationI ()2011af ). we recover a = 0.96. 
Since we find no evidence for non-gaussian scatter in- 
trinsic scatter in the Ysz — Yx relation, and furthermore, 
the observed scatter is dominated by measurement uncer- 
tainties, we believe that using the maximum likelihood 
fit is better justified. 

The final data set th at we compare against is that of 
lAndersson eralj poH . who find \nYsz/Yx = -0.20 ± 
0.1 using Chandra X-ray data and SPT SZ data with a 
median cluster redshift z « 0.7. This value is essentially 
in perfect agreement with ours, and is also consistent 
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Fig. 5. — The Ygz/Yx ratio, as observed using a variety of data. 
The horizontal band is our X-ray predictions model from section 
[21 the width of which (±0.05) reflect our ignorance of clumping 
corrections. The "PI. All" da ta point shows our re-analysis of the 
IPlanck Collaboration! Il2011al) data. "PI. Sub" uses this same data, 
but restricts the cluster sample to that e mployed in our Chandr a 
analysis. The Andersson point comes from lAndersson et al.l 112010 '). 



with that of iPlanck CoUaborationl ()2011a|) . The uncer- 
tainty in their measureme nt is significantly larger than 
the one reported here or in IPlanck Collaboration! (|2011a[ ) 
both because of a smaller cluster sample, and because 
the X-ray data is also shallower in terms of number of 
X-ray photons. We conclude that from the comparison 
of Planck -\- Chandra and SPT-I- C/iandra samples, there 
is no evidence for evolution in the Ysz — Yx relation. It 
will be interesting to make the same comparison in the 
future as both projects acquire more galaxy clusters with 
high quality SZ and X-ray data. 

For completeness, we have also estimated the am- 
plitude of the Ysz-Yx relation setting a = 1 for 
both our data, and the va rious sub-samples of the 
IPlanck Collaboration! (|2011al ) data we considered here. 
The likelihood for each of these fits is shown in Figure 
[SI In addition, we also computed the best fit amplitude 
for the low redshift (z < 0.12) and intermediate redshift 
("0.12 < z < 0.3) sub-samples from IPlanck Collaboration! 
(|2011a! l. Our results are summarized in Figure [5] 

6.2. Comparison to X-ray Expectations 

Our empirical calibration of the Ysz — Yx rela- 
tion for z < 0.1 Planck clusters is somewhat lower 
than the expectation from the observed cluster tem- 
perature profiles for both Chandra and XMM mea- 
surements, Alii{Ysz/Yx)xMM = -0.09 ± 0.035 and 
AlniYsz/Yx) Chandra = -0.127 ± 0.024. Some of 
this offset can be explained by small-level nonunifor- 
mities of the intracluster gas (gas clumping factor Q 
in eq. [5]). As per our earli er discussion, we set Q = 
1.05±0.05 (jNagai et al.ll2007[ ). which brings the predicted 
Ysz/Yx ratio within a few percent of the observed value, 
Alii{Ysz/Yx)xMM = 0.09 for the fuU XMM^Planck 
cluster sample, Ahv{Ysz IYx)xmm = -0.042 ± 0.035 
for the Chandra sub-sample of galaxy clusters, and 
A\n{Ysz I Yx) Chandra = -0.078 ± 0.024 as estimated 
from Chandra. Note that the estimated systematic un- 
certainty in the clumping correction is ±0.05, so all val- 



ues arc in reasonable agreement with expectations, even 
though some are in strong statistical tension with one 
another. 

Our results fit within the general picture that has de- 
veloped in the recent past. A few years back, a vari- 
ety of works argued that the SZ signal of galaxy clus- 
ters was inconsistent with predictions fro m X-ray data 
(|Lieu et al.! [2006t iBielbv fc ShmAsI [2007h . These au- 
thors fit isothermal /3 models to ROSAT data in or- 
der to predict the SZ-signal of galaxy clusters as ob- 
served by WMAP, finding that the observed SZ sig- 
nal was significantly lower than that their predictions. 
There is now ample evidence that (3 models are in- 
adequate fits to the intra-cl uster gas of galaxy clus- 
ters fe.g lVikhlinin et al.! [20061: iCroston et a l.il2008j). and 



that clusters are not i s othermal (e.g. I Vikhlini n et al 
2006| [Pratt et all 12007!: iSun et al.[ [2009: Morctt i et al 
20111 ). [Hallman et al.! (|2007[ ) and iMroczkowski et al 



explicitly demonstrated that isothermal model 
fits to X-ray data leads one to over-predict the SZ sig- 
nal of galaxy clusters. Using more realistic density 
and temperature profiles, one finds good agreeme nt be- 
tween SZ and X-ray observa b les ( A trio-Barandcla et aLl 
2008!: [Komatsu et al! [20081 iMroczkowsk et al.[ [20091: 



Diego &: Partridge! [201 Ol:"IZwart et al. ! [201 It [Melin et al.! 
2011!: [Planck Collaboration[[2011a![cD . The results sum- 



marize above fit well within this new wave of works. 
7. SUMMARY AND CONCLUSIONS 



Usin g SZ data from Planck (jPlanc 



c Collaboration 



2011a!) and X-ray data from Chandra ( Vikhlinin et al 
20091 ). we have estimated the Ysz-Yx scaling relation 



Assuming a power-law of the form in equation ll2i we find 
\n{Ysz/Yx) = a = -0.202 ± 0.024, a = 0.916 ± 0.032, 
and asz\x = 0.083 ± 0.021. The scatter constraint is 
prior-driven, and the variance in the relation is consis- 
tent with zero intrinsic scatter. The slope a differs from 
unity at the sa 2.3(7 level, which is interesting but not 
yet significant. We find no significant difference in the 
scaling relations when comparing relaxed and unrclaxed 
systems. 

Our mean Ysz/Yx is significantly lower than that re- 
ported from the XMM and Planck analysis of the full 
sample of 62 clusters (a = —0.032 ± 0.028). However, 
when we repeat the Planck and XMM analysis restrict- 
ing ourselves to the same 28 systems that compose the 
Chandra subsample, we find results that are in signifi- 
cantly better agreement (a = — 0.166 ± 0.036), demon- 
strating the original tension is driven by the remaining 
34 (mostly z > 0.1) clusters in the full Planck and XMM 
cluster sample. 

Compared with the XMM-Newton analysis, we find 
a significantly lower scatter in the Ysz — Yx relation, 
<^sz\x ~ 0.083 ±0.021 for the Chandra- Planck data, ver- 
sus a^z\x = 0.182 ±0.025 for the full XMM-Planck sam- 
ple. The scatter for the full XMM-Planck sample is 
at least partly driven by the systematic differences be- 
tween the Chandra sub-sample of galaxy clusters, and 
the remaining systems. The XMM-Planck scatter for 
the Chandra subsample is lower asz\x = 0.132 ± 0.036 
and consistent with the scatter estimated from Chandra 
data. Finally, we note the Chandra data is fully con- 
sistent with no intrinsic scatter in the Ysz~Yx relation: 
our central values are prior-driven. 
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Both the Chandra and XMM values of Ysz/Yx for 
z < 0.1 are significantly lower than a — —0.075 ex- 
pected from the observed cluster temperature profiles 
along (§ EJ- We believe that a significant fraction of 
that difference can be attributed to some level of gas 
clumping, with Q 1.05 including substructure masking 
(|Nagai et al.ll2007[) . Any remaning difference is within 
the range of current absolute accuracy of the X-ray tem- 
perature measurements and systematic uncertainty in 
the clumping factor correction. 

In summary, we find that the observed Yx — Ysz re- 
lation from the joint Planck and Chandra analysis is in 
good agreement with the theoretical expectations — it 
shows very low intrinsic scatter, a slope close to unity, 
and the amplitude is consistent with physically plausible 
levels of gas clumping. 

In conjunction with the work in lVikhlinin et all (j2009f ). 



our results can be used to self-consistently estimate the 
multi-variate scaling relation of galaxy clusters. In an 
upcoming work, we use these results to study the Lx-M 
and Ysz-M scaling relations of galaxy clusters as de- 
termined from X-ray data, and compare these to those 
derived fror a the optically select ed maxBCG galaxy clus- 
ter catalog (jKoester et al.ll2007[ ). 
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APPENDIX 
CLUSTER DATA 
We summarize the data employed in our analysis in table [TJ 
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Fig. 6. — Best-fit parameters for the Ygz^^x scafing relation assu ming a slope of unity. Black c urves are obtained using Chandra data 
and the cluster sample in Table [T] Red and b lue curves are using thc lPlanck Collaboration! I l2011al ') data. The red curves show the results 
when we restrict our analysis of the IPlanck Collaboration (2011a) data to the same galaxy clusters that constitute our Chandra sample. 
The blue curves include all galaxy clusters in iPlanck Collaboration. (,2011al V All confidence regions are 68% and 95% likelihood contours 
respectively. 
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TABLE 1 
Cluster Data 



Name 


Pll Rcdshift 


Rcdshift 


A 


CYx {Chandra) 


CYx (XMM) 


Merger? 


ASS 





.052 





.0557 


5.05 ±0.54 


6.19 ±0.16 


5.37 ± 0.22 


— 


A119 


0. 


.044 


0. 


.0445 


2.60 ±0.29 


3.55 ±0.13 


3.42 ± 0.16 


/ 


A401 


0. 


.075 


0. 


.0743 


7.71 ±0.74 


10.86 ±0.47 


10.42 ±0.75 


— 


A3112 


0. 


.070 


0. 


.0759 


1.98 ±0.33 


3.38 ±0.21 


2.82 ± 0.11 


— 


A3158 


0. 


.060 


0. 


.0583 


3.14 ±0.27 


3.10 ±0.07 


3.73 ± 0.15 


— 


A478 


0. 


.088 


0. 


.0881 


8.71 ±0.76 


10.17 ±0.58 


9.59 ± 0.39 


— 


A3266 


0. 


.059 


0. 


.0602 


8.84 ±0.69 


11.89 ±0.29 


10.07 ±0.36 


/ 


A3376 


0. 


.045 


0. 


.0455 


0.97 ±0.19 


1.64 ±0.24 


1.34 ±0.06 


/ 


A1413 


0. 


.143 


0. 


.1429 


6.99 ±0.76 


7.83 ± 0.48 


7.58 ± 0.12 


— 


ZwC11215 


0. 


.077 


0. 


.0767 


4.32 ±0.66 


5.54 ±0.21 


5.72 ± 0.30 


— 


A3528S 


0. 


.053 


0. 


.0574 


2.41 ±0.33 


2.00 ±0.15 


1.62 ± 0.10 


— 


A1644 


0. 


.047 


0. 


.0475 


2.41 ±0.39 


3.08 ±0.30 


2.80 ±0.13 




A1650 


u. 


.Uo4 


U. 


nQOQ 

.Uozo 


4.00 ±0.55 


3.62 ±0.15 


3.67 ±0.08 


— 


A1651 


0. 


.084 


0. 


.0853 


3.50 ±0.58 


5.32 ±0.29 


4.12 ± 0.12 


— 


A1689 


0. 


.183 


0. 


.1828 


12.93 ± 1.42 


13.42 ±0.71 


12.41 ±0.22 


— 


A3558 


0, 


.047 


0. 


.0469 


3.95 ±0.47 


4.04 ±0.36 


4.50 ±0.18 


/ 


A1795 


0, 


.062 


0. 


.0622 


4.37 ±0.38 


5.34 ±0.18 


6.78 ± 0.28 




A1914 


0, 


.171 


0. 


.1712 


9.47 ±0.85 


13.29 ±0.68 


12.43 ±0.31 




A2034 


0, 


.151 


0. 


.1130 


4.27 ±0.58 


5.11 ±0.15 


6.01 ± 0.14 




A2029 


0, 


.078 


0. 


.0779 


7.65 ±0.66 


11.48 ±0.53 


12.13 ±0.65 




A2065 


0, 


.072 


0. 


.0723 


3.72 ± 0.48 


4.35 ±0.16 


4.52 ± 0.23 


/ 


A2163 


0, 


.203 


0. 


.2030 


43.01 ± 1.99 


60.93 ±3.21 


59.75 ± 2.14 


/ 


A2204 


0, 


.152 


0. 


.1511 


10.39 ± 0.94 


12.96 ± 1.00 


11.88 ±0.39 




A2256 


0, 


.058 


0. 


.0581 


6.73 ±0.38 


8.94 ± 0.43 


7.02 ± 0.33 


/ 


A2255 


0, 


.081 


0. 


.0800 


4.81 ±0.37 


4.97 ±0.14 


4.81 ± 0.15 




RX J1720 


0, 


.164 


0, 


.1593 


5.68 ±0.72 


6.45 ±0.38 


5.69 ± 0.14 




A2390 


0, 


.231 


0. 


.2302 


15.61 ± 1.22 


19.49 ± 1.56 


19.26 ±0.58 




A3921 


0, 


.094 


0. 


.0941 


3.41 ±0.28 


3.75 ±0.15 


3.10 ±0.08 





Note. — The units for D\Ysz and CYx are 10"^ Mpc^. The Yx values from XMM are drawn 
from | Piancl<: Collaboration! Il2011al) . In addition, both the Ygz and Yx data of lPlanclt Collaboration! 
I l2011al ) are corrected to our fiducial cosmology (h = 0.72) and to our cluster reds hifts. The Pll 
redshift column quotes the redshifts as they appear in IPIanck Collaboration! Il2011al 'l. 



